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■ A quantum generalization of the semiclassical theory of Gutzwiller is given. 
The new formulation leads to systematic orbit-by-orbit inclusion of higher Ti 

\ contributions to the spectral determinant. We apply the theory to billiard 

0^ ' systems, and compare the periodic orbit quantization including the first h~ 

■ contribution to the exact quantum mechanical results. 
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The Gutzwiller trace formula jl| is the most compact formulation of the semiclassical 
quantization of classically chaotic systems. It expresses the trace of the Green function in 
terms of classical periodic orbits of the system. In recent years it has been demonstrated 
]l||J] on many classically chaotic systems that this is indeed a very good approximation. 
Nevertheless, it is only the leading order formula in h and systematic inclusion of higher 
order h contributions is desirable. In this letter we present a new method, which simplifies 
the computation of these contributions and gives some insight in their structure. Our key 
result is the generalized Selberg product ([14]) which expresses the global eigenvalues in terms 
of local eigenspectra computed on neighborhoods of individual periodic orbits. 

The trace of the energy domain Green function is 

TrG(x,x',E)= [ dxG(x,x,E)=J2ir^r, (1) 

where E n are the eigenenergies of the system. The trace is a logarithmic derivative of the 
spectral determinant 

TrG(x,x',E) = -^logA(E), (2) 

where A(E) = det(.E — H) = Yl n (E — E n ) is the formal secular determinant of the Hamilton 
operator H. The Green function can be written || as a sum over classical trajectories 
connecting x and x' with energy E 

G(x, x', E) = Y J A(x, x', E)e lS ^ x '' E)/h , 

cl 

where S(x, x', E) is the classical action of the trajectory, the amplitude A(x, x', E) can be 
determined by the multiple scattering expansion. For h — > the trace of the Green function 
can be evaluated with the saddle point method. The saddles of the function S(x, x, E) are 
the classical primitive periodic orbits p with energy E, and their repetitions r. Consider the 
trace of the Green function restricted to a small tube Q p enclosing the periodic orbit p: 

Tr G P {E) = drdq 1 dq 2 ...dq d _ 1 G(r, q, r, q, E), 

where r is a local coordinate along the orbit and q^i = 1, ...,d— 1 are coordinates or- 
thogonal to the orbit. We shall expand the action S Pir (r, q, r', q') and the amplitude 
Ap )r (r, q, t', q') around the periodic orbit and repetitions in Taylor series and TrG p (E) 
can be written as a trace of the locally expanded Green's function G p (t, q, r', q', E) = 
J2 P ,r A Pt r(T, q, t', q', E)e tSv ' r{ - T,cl ' T ' q ' E ^ h . The trace of the full Green function can be expressed 
with these local traces 

Tr G(E) = TiG (E) + ]TTrG p (£), (3) 

v 

where TrGo(E) is the Weyl or Thomas- Fermi term || coming from zero length trajectories. 
For h — > the values of the integrals are independent of the sizes and overlaps of the fl p 
tubes. As we shall be expanding around the classical orbits, our results are valid only in the 
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sense of asymptotic series; finite % computations require care, especially for non-hyperbolic or 
non-elliptic situations where periodic orbits might accumulate. In applications the number 
of periodic orbits used should depend on % and on the order of truncation of the Taylor 
expansion of the action and amplitude. Berry and Keating || have proposed a method 
for dealing with this problem; here we assume that our computations (here only the first h 
contribution will be evaluated) are within the scope of the validity of saddle point expansions. 

The standard way of evaluating the above integrals is the Feynman diagram perturbative 
exapansion about the saddle point, implemented in ref. jn|. However, this is not the only 



way to compute the trace of G p ; here we offer an alternative, based on computing the 
local eigenspectrum. Since G p is the locally expanded form of the Green function of the 
Schrodinger equation, we can assume that there exist a suitable local expansion of the 
Schrodinger equation whose Green function is G p . If we can solve the local Schrodinger 
problem on Q p , then the trace of G p can be recovered as a logarithmic derivative of the 
spectral determinant A p (E) of the local eigenvalue problem 

TiG p {E) = ^\ogA p {E). (4) 

To define the local Schrodinger problem we introduce a representation of the Schrodinger 
equation which is closely related to classical dynamics and facilitates development of a 
semiclassical expansion around a given classical orbit 0. We take an arbitrary solution of 
the Hamilton- Jacobi equation 

d t S(x, t) + \(VS{x, t)f + U(x) = 0, (5) 

and transform the wave function ip(x,t) of the Schrodinger equation to (p(x,t) = 
tp(x,t)e- iS( - x > t)/h . The new wave function satisfies the Maslov-Fjedoriuk transport equation 

dtcp + V( V VS) - ±<pAS - yAy? = 0. (6) 
The quantum mechanical propagator in this representation is 

(x', t'\x, t) = e £(s(^)-s(-'/))£( x ' 5 t >\ Xj t ) } (7) 

where C(x' ,t'\x,t) is the propagator of the transport equation @. By the analogy between 
(§) and the Fokker-Planck equation we can write this propagator as a sum over stochastic 
trajectories fluctuating around classical orbits 

C(x',t'\x,t) = (eiSf™WM^6(J-W,t,x))) 

where £(t',t,x) is the solution of the complex Langevin equation 

§ = V£(£, + ^(0 (8) 

with initial condition £(t) = x, and rj(t') is Gaussian white noise with (f]i(t)) = and 
{Vi(t)rij(t')) = 5ij5(t — t') . The symbol () denotes averaging over all noise configurations. 
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This representation based on an arbitrary solution of the Hamilton- Jacobi equation has 
been proposed in ref. || as an alternative to the Feynman path sum [fL0|1 . As H — > the 
fluctuations are suppressed, (§) becomes a deterministic equation and £(t',t,x) becomes the 
classical trajectory with initial coordinate x and momentum p(t) = VS(x,t). We develop 
a multiple scattering expansion of the propagator by handling the diffusive term — in 
(|) perturbatively 

oo 

C(x',t'\x,t) = C i0) (x',t'\x,t) + ]T(^/2)™C (n) (:r'^'lM), 

n=l 

where 

£<°V. t'\x, t) = ^W'\t, x ),n<*" 5{x > _ ^ tj x)) 

and £(t',t,x) is the deterministic trajectory. To study the propagator (0) for x, x' G Q p 
we can choose S p (x,t) such that the momentum p(t) = VS p (x p (t),t) coincides with the 
momentum of the periodic orbit p p (t) and the velocity field VS p (x, t) is independent from t. 
Such a local velocity field can be constructed order by order by substituting the the Taylor 
expansion 

oo 

S p (x,t) = Y / Sn(t)(x - x p (t)) n /n\, (9) 

n 

into the Hamilton- Jacobi equation (^). The symbol n = (n 1; n2, n^) denotes the multiple 
index in d dimensions, n! = Y\f =1 ni\, and (x — x p (t)) n = Y\i = i(xi — x Pti (t)) ni . Due to the time 
independence of the velocity field, the expansion coefficients are periodic functions s n (t) = 
Sn(t + T p ), except the zero order term which grows with each period by So(t + T p ) — So(t) = 
Jq p (^Pp(t ! ) — U(x p (t'))^j dt' = Jq p L(t')dt' . This solution is stationary d t S p (x,t) = —E due 
to the time independent velocity field. By order-by-order construction in each step we 
get a formally periodic velocity field, therefore it is advantageous to study the eigenvalue 
problem using the one-period operator ( 0) p , where p indicates the computation of 

the propagator with the above procedure using S p (x,t). If ip p (x,t) is a stationary state of 
the local problem, then it is also an eigenfunction of the one-period propagator 

e- iETp/h tfj p (x,0) = J dx'(x,T p \x',0} p ip p (x',0). (10) 

Using the fact that S p (x,T p ) — S p (x, 0) = J P dt L(t') and that the action of the periodic 
orbit is S p (E) = § p p dq = J P L(t')dt' + ET p we get for the transformed wave function 

<p p (x, 0) = e iSp{E)/h J dx'Cpix, T p \x' } 0)<p p (x', 0). (11) 

If the eigenvalues e Xp ' m ^ of C p are known 

e Vm(i? V P ,m(x, 0) = J dx'C p {x, T p \x', 0)ip p , m (x', 0), (12) 

where m is the multiple index labeling the eigenstate, the local spectral determinant is given 
by 
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A P (E) = - e x ^ E) e^ E) ). (13) 

m 

Now we can reexpress the spectral determinant of the full problem through an d 



A(E) = e w ^ I] IK 1 " e x ^ei 3 ^), 



(14) 



p m 



where W(E) = J E dE' Tr Gq(E') is the Weyl term contribution. This Selberg-product type 
expression for the spectral determinant is our main result; it is the full quantum general- 



ization of the semiclassical Gutzwiller-Voros product formula [p^] . A similar decomposition 
has been found for quantum baker maps in ref. The functions 



(15) 



are the generalizations of the Ruelle dynamical zeta functions flTa] . The corresponding 
generalization of the semiclassical Gutzwiller trace formula to the full quantum mechanical 
trace follows from (0); it can be writen as a sum over the prime cycle repetition index r: 



TyG(E) = TtG (E) 



1 

ih 



E 

p,r,m 



T P (E) - ih 



dXp, m (E) 
dE 



(E) e j-S P (E) 



(16) 



In practice the eigenvalues of ([12]) can be computed the following way. We substitute 
the perturbation expansion (p Ptm {x, t) = J2iLo (y ) V^mO^i t) i n t ne equation (|) and we get 
an iterative scheme starting with the semiclassical solution y^ ^: 



oHV + MV^Sp + l^ ^p 



(17) 



,(<+l)AO 

The eigenvalue can also be expanded in powers of ih/2: 



m" 



\,m{E) = 



ih 



C (i) 

p,m 



(18) 



1=0 



The eigenvalue equation (|T2| ) in ft expanded form reads as follows 



(0) 



,T p )=exp(cWj^(x l 0) l 

, T p ) = exp(C(°) 1 ) {^(x, 0) + Cfti^s, 0)} 



(19) 



and so on [|Tl|] . This approach is very similar to the one used in ref. |T2j . 

We expand the functions (pf^ix, t) as in (^), in Taylor series around the periodic or- 
bit, and to solve recursively the equations (|T8|). Only a couple of coefficients should be 
computed to derive the first corrections; the detailed calcualtion will be published else- 
where |Tl||. The zeroth order term yields the Gutzwiller semiclassical weight = 
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—iiriSp — Yh=i ( m « + |) u p,ii where u Pji = log |A Pji | are the Lyapunov exponents of p, \A Pji the 
eigenvalues of the monodromy matrix J p , and v p is the Maslov index of the periodic orbit. 
The first correction is given by the integral C^ m = Jq p dt Aip^ m (t) / (p^ m (t) . 

The semiclassical trace formula can be recovered by dropping the subleading term 
-ihjgXp^E) and using the semiclassical eigenvalue e A p° m ^) = = e ~ lUp7T ~^ ( - mt+1/2 * >Up - 1 . 
Summation for the indexes m; yields the semiclassical amplitude 

rni p -irupw 

Y e x v ^ E)r = - (20) 

|det(i-J;)|V 2 - W 

When the theory is applied to billiard systems, the wave function should fulfill the 
Dirichlet boundary condition on hard walls. The wave function determined from (|6]) behaves 
discontinuously when the trajectory x p {t) hits the wall. For simplicity we consider a two 
dimensional billiard system here. The wave function on the wall at the time t-o immediately 
before the bounce is given by 

1> in (x, y{x), t) = <p(x, y(x),t_ )e iS ^^-°y h , (21) 

where y(x) = Y 2 x 2 /2\ + y 3 x 3 /3! + Y±x A /<i\ + ... is the parameterization of the wall around 
the point of reflection. The wave function on the wall at the time t + o immediately after the 
bounce is given by 

iM*, y(x),t) = <p{x, y(x),t +0 )e lS ^^ n . (22) 

The sum of these wave functions should vanish on the hard wall. This implies that the 
incoming and the outgoing amplitudes and the phases are related as S(x,y(x),t_ ) = 
S(x,y(x),t +0 ) and (p(x, y(x), £_o) = — f(x, y(x), t+o). The minus sign can be interpreted 
as the Maslov phase associated with the hard wall reflection. 

In order to gauge the improvement due to the inclusion of the quantum corrections we 
have developed a numerical code [|1(J for calculating the first correction Co 1 ^ for general 



two-dimensional billiard systems. The first correction depends only on the basic periodic 
orbit data such as the lengths of the free flights between bounces, the incidence angles and 
the first three Taylor expansion coefficients 13, 14 of the wall in the point of incidence. 
To check that our new local method gives the same result as the direct calculation of the 
Feynman integral, we computed the first h correction C p ^ for the periodic orbits of the three 
disk scattering system ]T7|JT3. We have found agreement to the five decimal digits listed 
for the quantum corrections calculated in ref. |T9 (our method generates these numbers to 



any desired precision). The m 7^ coefficients cannot be compared to ref. [19], since the 
m dependence was not realized there due to the lack of general formulas (14) and flT5|). 
However, the m dependence can be checked on the 2 disk scattering system | 20| . For the 
standard example [p^-pO[], with the distance of the centers (R) taken to be 6 times the disk 
radius (a), we obtain = (-0.625m 3 - 0.3125m 2 + 1.4375m + 0.625)/v / 2~E. For m = 



and 1 this has been confirmed by A. Wirzba pl| , who was able to extract these terms 
from his exact quantum calculation. Furthermore, our calculation has been verified for the 
confocal hyperbola system at many different parameters p2|. Our method makes it possible 



to utilize the symmetry reduction and to include quantum coreections into the fundamental 



domain cycle expansion calculation of ref. [18]. We have computed the first correction 
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to the leading 226 primitive periodic orbits with 10 or less bounces in the fundamental 
domain. Table I demonstrates that the error of the corrected calculation vs. the error of the 
semiclassical calculation decreases with the wave number. Besides increasing the accuracy, 
a fast convergence up to six decimal digits can be observed, in contrast to the three decimal 
digits obtained in the full domain calculation [19| . 
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manuscript. We thank P. Cvitanovic the conceptual reedition of our manuscript. G. V. 
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TABLES 



TABLE I. Real part of the resonances (Re k) of the three disk scattering system at disk 
separation 6:1. Semiclassical and first corrected cycle expansion versus exact quantum calculation 
and the semiclassical error Ssc = |L- e (^QM) — Re(&sc)| divided by the error of the first correction 
Scorr- The magnitude of the error in the imaginary part of the resonances remains unchanged. 
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u bC 1 u Corr 


fl 6Q7995 


D 75831 3 


fl 5851 5D 


n 53 


2.239601 


2.274278 


2.222930 


2.08 


3.762686 


3.787876 


3.756594 


4.13 


5.275666 


5.296067 


5.272627 


6.71 


6.776066 


6.793636 


6.774061 


8.76 


30.24130 


30.24555 


30.24125 


92.3 


31.72739 


31.73148 


31.72734 


83.8 


32.30110 


32.30391 


32.30095 


20.0 


33.21053 


33.21446 


33.21048 


79.4 


33.85222 


33.85493 


33.85211 


25.2 


34.69157 


34.69534 


34.69152 


77.0 
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